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The most dramatic variation across the superconducting transition temperature 
(Tc) in spectral function is the development of peak- dip -hump (PDH) structure near 
(7r,0)n!|. This has been observed by angle-resolved photoemission spectroscopy 
(ARPES) in many high temperature superconductors (HTSC's), and is thought to 
carry critical information about the superconducting transitionpi. Particularly, on 
Bi2Sr2CaCu208+5 (Bi2212), we found the superconducting peak (SCP) intensity to be 
closely related to the condensate fraction of the system This was achieved by a 
spectral analysis based on the SCP ratio (SPR) accounting for the phenomenologically- 
extracted coherent component in (tt, 0) ARPES lineshape. The phase-coherent like 
behavior of SPR (increases with doping increases) in the underdoped regime was later 
confirmed by H. Ding et al. based on a similar analysis of the lineshape by means of the 
so-called renormalized quasiparticle weight (RQW), a quantitive equivalent to SPRjlj. 
However, as a contrast, the decrease of the SPR upon doping increase across the optimal 
doping was not reproduced. This led the authors to think it was necessary to combine 
the additional effect by the gap opening to yield, rather than simply taking the SPR as, 
the long-sought supercoducting order parameter for the cuprates. 

This controversy remains until later another important finding was made, the 
existence of bilayer band splitting (BBS) in heavily overdoped (OD) Bi2212|Hj. This 
gives birth to a sharp anti-bonding band remnant (AB hump) even in the normal 
state at the SCP binding energy position tending to invalidate the SPR (RQW)-based 
phenominological background substraction: the SPR (RQW) obtained this way contains 
a substantial spectral weight contribution from the AB hump, which should not be 
counted in in order to obtain a true ratio reflecting the SCP only. The increasing 
overestimation of the true peak ratio prevents a similar analysis from applying to 
the more OD regime where the AB band becomes sharper. In Ref. [6], a lineshape 
modelling based on a physically sensible description of the two quasi-particle bands [7] 
proved its eligibilty to reproduce the complex PDH lineshape under a considerable BBS 
complication and yielded a reduced SPR strengthening its correlation with the superfiuid 
density in the very OD regime. 

In this paper, we revisit the investigation by presenting more results reflecting the 
robustness of the conclusion, which include the insensitiveness to a physical variation in 
the spectral function, to a different self-normalization scheme adopted. 

For clarity, we recap below the fitting function used for the lineshape modelling in 
Refs. inilHl- Note that only possible variations in some expressions are pointed out and 
the denominations of symbols completely follow those in Refs. [HI El otherwise. 
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Figure 1. (color) Results of Scheme II of the global convolve fit on the T-dependent 
(a) (see also [f] for an enlargement of the T=30 K curves) and /iz^-dependent ARPES 
spectra at (tt, 0) of (Pb)-OD65 in the superconducting state at T=10 K (b) and normal 
state at T=80 K (c). The AB hump lies lower in energy than the BB hump. The width 
difference of the SCP is due to resolution variations from 10 to 18 meV at different 
hv's. The fitted /iz/-dependent matrix elements (d) and T-dependent energy positions 
of the maxima (e) of various spectral features are shown reminiscent of the similar 
results in Ref . [S] ■ Inset of (e) illustrates the constancy of the AB peak width upon T 
variation in the superconducting state. 
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where the imaginary part of the self-energy is defined either by 



S"(^,T) = v'(a^)2 + (/5T)2 (1) 
, as in Ref. P, or by 

E"{uj, T) = I ^{auy + {f3Ty + CI (2) 

[S], where hfetime broadening associated with the impurity level of samples is reflected 
by ( and that with doping level of by a and (5, which are photon energy {hu) and 
temperature (T) independent. r(T) is a T-dependent linewidth (FWHM) for the SCP. 
Note that, as a sum rule for an overall spectral weight conservation upon T variation, 
we adopt either 

/oo 
Ah(p){uj)duj = 1 (3) 
-oo 

Refs. P IHI or 

/■oo 

f{uj)Af,^p){uj,T)duj = n{T) (4) 



PP, where n{T) is the occupation number at (vr, 0), which is nearly independent of 
and thus assumed to be a constant here. As a result, spectral weight transfer (SWT) 
(= 1 — Ca(b), < Ca(b)(T) < 1) is well defined for the description of the SCP formation 
disregarding the kind of dependence under investigation p. In any case, a realistic 
SCG or pseudogap opening extracted with the leading edge gap method on the same 
sample introduces an effective T-dependent energy shift to the bilayer-split bands in the 
superconducting state while the energy position of the SCP is subjected to the fit. 

In order to retrieve reliable quantitative information from a multi-dimensional 
fit to a set of interconnected EDC's, we developed a powerful fitting routine under 
WaveMetrics Igor Pro 4.0 [TUj, which has been applied in our previous workp and 
detailed in Ref. jH]. Physically, the goal is to achieve a systematic understanding 
of both the hi'- and the T-dependent behavior of the lineshape where many physical 
constraints are involved and strictly followed. For example, in the the /izz-dependent set, 
all functions of T are shared variants to all EDC's at the same T while the functions of 
hu (matrix elements) are locally specified, and in the T-dependent set, vice versa. The 
renormalized band energy position e is always a global fitting parameter for the EDC's 
of samples at the same doping level. Technically, reminiscent of the great contrast in the 
overall quality (or details) of the fit between the one shown in Fig. 2 of Ref. [7] and ours, 
in Refs. P P, though on data sets generally twice larger in size, the inherently local 
fit scheme [71 ITT] with manual adjustments on the globally defined fitting parameters is 
unlikely to give a reliable result in a truly global sense. In contrast, our fully-automatic 
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Figure 2. (color) Results of the three schemes of the global non-analytical fit. The 
fitted /ii/-dependent matrix element hump ratios, (filled markers), and peak- 

hump ratios, ^^+'Mbh)/2 (''^P^^ markers) are shown in (a), (b) shows the fitted 
T-dependent energy positions of the AB peak (empty markers) and AB hump maxima 
(filled markers), (c) The T-dependent SWT's and an interpolation on the averaged 
dependence according to the function SWT=SWTo[l - (T/TcY] (see Ref. US]). The 
fitted value, C = 2.56, agrees with the T-dependent /xSR result on Tl2Ba2CuOQ+s 
at the same doping level (OD66 T12201)(not shown). The contrast between the 
phenomenological SPR and SWT-based SPR's extracted from the same set of data 
(taken with He-I) is given in the inset, (d) The doping-dependence of SPR {hv~22.7 
eV, T = 10 K) supplemented by the global fit results on Pb-OD65. 

fitting routine can yield a robust, physically-constrained global fit for a large set of 
EDO's regardless of the initial values input. 

The global non-analytical fit implemented in Ref. [B| is based on Eq. ^ and El 
(Scheme 1). In the following, we conducted another two schemes, the one with Eq. |2l 
andini (Scheme 11) and the other with Eq. |21andE]to complement the results of analysis 
based on the global non-analytical fit. Our argument about the qualification of the SPR 
as the order parameter for the HTSO's is invulnerable to the physical variation of the 
fitting scheme we choose. This checkup is hardly practical without a fully-automatic 
multi-level global non-analytical fit routine which can be easily oriented for any non- 
analytical global fitting task given. 

In Fig. ^ we present the results of Scheme 11, reminiscent of Fig. 2-3 in Ref. [H]. 
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The fitting quality is exemplified in Fig. WIS)- Keep in mind that this quality is common 
to the fits of the 35 EDO's in contrast to a local best fit on a few selected curves at 
the expense of a global agreement on a much larger assembly within the usual manual 
fitting framework[7]. In Scheme II, when we further release the constrain on Ta to let 
it fittable, we find its T-dependence is negligible as obtained by our phenomenological 
fitting. This argues against the Fermi liquid approach that attributes the disappearance 
of the sharp peak above Tc to broadening caused by phase fluctuations [T^ IT^ . 

Note that the addition of an impurity term (whose fitted value is —0.0058 ±0.0004 
for Pb-OD65 and —0.0067 ± 0.0004 for OD65) in the imaginary self-energy doesn't 
result in any qualitative change in the /iz/-dependence of matrix elements as well as 
the T-dependence of the energy positions of spectral features [Fig. [Ud)(e)]. However, 
quantitative changes is discernible when we plot the results of the three schemes together 
as in Fig. El Though the modification of the normalization equation is not expected 
to alter much the energy positions of features whose maximums are more or less self- 
pinned [Fig. Efb)] and the anti-phase-like /iz^-dependence of the hump matrix elements 
as predicted by theory [121 En], the matrix element ratio of the two hump bands is, 
in principle, a touchstone for the justification of the physical normalization method. The 
average of matrix element hump ratio over a large hu window is supposed to approach 
unity as expected from a pair of bands of the same atomic character split by the c-axis 
bilayer coupling. Though the hu window concerned here is not large enough. Scheme I 
and II give a ratio of 0.91 ± 0.01 and 1.21 ± 0.04, respectively, quite approaching unity, 
while Scheme III yields 1.55 ± 0.02. 

Interestingly, compared with the hump ratio, the matrix element ratio between the 
peak and the hump, defined by ^j^^ ^^2/6^)72 ' exhibits a constancy over hu. Reminiscent 
of the absence of SCP finding report in spectra of mono-layer cuprate Bi2Sr2Cu06+5, 
this probably imply the AB peak has a hybrid origin of AB and BB humps at the same 
time. It is argued by the magnetic resonance (vr, vr, vr) mode scenario that the c-axis 
momentum transfer by the mode scattering occurs between AB and BB which leads 
to the formation of the PDH lineshape^Zj- However, in our picture the SCP is a new 
quasi-particle developed from its normal state counterpart of hump, just the opposite 
to the magnetic mode scenario where the dip (hump) is argued to contain the incoherent 
spectral weight subjected to the scattering of the quasi-particle peak. This critical 
question concerning the SC origin is still unaddressed. 

As told by Fig. Efc), the SWT saturates at low temperature to be (0.26 ± 0.03), a 
value far below 1, which is required by the magnetic mode scenario in which the collective 
mode induces an on/off damping mechanism at the dip position driving the whole lower- 
lying AB component coherent with no AB hump remanence|18j. As mentioned before, 
the large remanence of AB hump proportion at the experimentally observed SCP energy 
position in the superconducting state also questions the validation of the SPR (RQW)- 
based phenomenological background substractionO |3] in the very OD regime. This 
can be illustrated by the contrast in the inset of Fig. Wic) with the phenomenological 
SPR featured by an overall overestimation in value and a non-zero normal state baseline 
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followed by a rounded onset spreading over a large T region across T^. 

The collapses of the T-dependent SWT (SPR) onto the superfiuid density 
curve [Fig. Efc)] and of the SWT-based SPR's by the three schemes onto the Tc 
guideline strength [Fig. Efd)] strengthen again the qualification of SPR (SWT) as a 
superconducting order parameter for Bi2212. However, the component-resolved matrix- 
element-free merit of SWT makes it an ideal quantity for the quantitative spectral 
analysis of the complex lineshape at (7r,0). A systematic investigation of both the 
doping and T dependence of SWT based on an indispensable global fit facility is needed 
for a comprehensive understanding of its relation to the possible mechanism for high- 
temperature superconductivity. 
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